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On the Lamperti stable processes 



M.E. Caballerc0 J.C. Pardfland J.L. Perez^ 



Abstract 

We consider a new family of IR'^-valued Levy processes that we call Lamperti stable. One of 
the advantages of this class is that the law of many related functionals can be computed explicitly. 
Here, in the one dimensional case we provide an explicit form for the characteristic exponent and 
other several useful properties of the class (for some particular cases see [6], [lOj . |24| and [28]). 
This family of processes shares many properties with the tempered stable and the layered stable 
processes, defined in Rosiiiski [31] and Houdre and Kawai [16] respectively, for instance their 
short and long time behavior. We also find a series representation which is used for sample path 
simulation, illustrated in the case d = 1. Finally we provide many examples, some of which appear 
in recent literature. 

Key words and phrases: Lamperti stable distributions and processes, stable processes. Levy 
processes. 

MSG 2000 subject classification: 60E07, 60G51, 60G52. 

1 Introduction. 

In recent years the interest in having more accurate models in various domains of applied probability 
has lead to an increasing attention paid to some special classes of Levy processes related to the stable 
law, for example: the tempered stable and the layered stable processes introduced by Rosinski in 
|31j and Houdre and Kawai in [16 , respectively. Both families of processes have nice structural and 
analytical properties, such as combining in short time the behavior of stable processes and in long 
time the behavior of a Brownian motion. They also have a series representation which may be used 
for sample paths simulation. 

Lamperti [25] and more recently, Caballero and Chaumont [6] studied four families of Levy pro- 
cesses which are related to the stable subordinator and some conditioned stable processes via the 
Lamperti representation of positive self-similar Markov processes. Those studies are the starting 
point of our work. Recall that positive self-similar Markov processes, {X,Px), x > 0, are strong 
Markov process with cadlag paths, which fulfill a scaling property, i.e. there exists a constant a > 
such that for any 6 > 0: 

The law of (6Xf,-at, t > 0) under is ¥bx- 

We shall refer to these processes as pssMp. According to Lamperti [5S], any pssMp up to its first 
hitting time of may be expressed as the exponential of a Levy process, time changed by the inverse 
of its exponential functional. Reciprocally, any Levy process ^ can be expressed as the logarithm of a 
time changed pssMp X. In this paper we refer to this Lamperti transformation as LTi and the details 
can be seen in [25]. 
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One of the examples treated by Lamperti in [2S] is the case when {X,Px) is a stable subordinator 
of index a G (0, 1) starting from x > 0. Lamperti in [2S], describes the characteristics of the associated 
Levy process s = (sj, t > 0) which is again a subordinator, with no drift and with Levy measure given 

by 

nidx) = — — rrdx, a; > 0. 

The three cases of pssMp studied in [B] are related to some conditioned stable processes. The first 
one is the stable Levy processes killed when it first exits from the positive half-line, here denoted by 
{X* ,Px)- The second class corresponds to that of stable processes conditioned to stay positive (see 
for instance [El [11]), denoted by {X\Px)- Finally, the third class of pssMp is that of stable processes 
conditioned to hit continuously, denoted by {X-^,Px)- The corresponding Levy processes under the 
LTi transformation are denoted by £,* ,£,^ and respectively. These three classes of Levy processes 
have no gaussian component and their Levy measure are of the type 

= (^ (e._l)a+i J^{^>o} + (1 _ ^{^<o} j dx, 

where c+, c_ are the constants of the Levy measure of the original stable process and 6 is a positive 
parameter. We recall that for ^* the constant b equals 1 and moreover it has finite lifetime if c_ > 
and, when c_ = 0, it has infinite lifetime and drifts to —oo. For the constant b is equal to ap+ 1, 
where p = Pq{Xi < 0). It has infinite lifetime and drifts to oo. Finally, for the processes the 
constant b is ap. It has infinite lifetime and drifts to — oo. We remark that such processes have linear 
coefficients in their respective characteristic exponent that we denote by a* , and . Such constants 
are computed explicitly in [^ in terms of a, p, c_ and c_|_. Actually it was proved recently in |10| . 
that the process corresponds to conditioned to drift to — oo (or equivalently is conditioned 
to drift to +00). 

Finally, motivated by self-similar continuous state branching processes with immigration, Patie 
has recently studied in [28 the family of Levy process with no positive jumps with Levy measure 

77*(dx) = J- r — rrdx, x < 0, 

where i9 > —a. 

These Levy processes have the advantage that the law of many functionals can be computed 
explicitly, for example: the first exit time from a finite interval or semi-finite interval, overshoots 
distributions and exponential functionals (see for instance Caballero and Chaumont [B], Chaumont 
et al. [10 , Kyprianou and Pardo [23] and Patie [H]). We also emphasize that in some cases the 
Wiener-Hopf factors can be computed and scale functions in the spectrally one sided case can be 
obtained. Since many tractable mathematical expressions can be computed, this class seems to be an 
useful tool for applications and rich enough to be of particular interest. 

In this work, we investigate a generalization of the Levy processes mentioned above and we will 
refer to them as Lamperti stable processes. We also study these processes in higher dimensions. We 
will see that this class has nice structural and analytical properties close to those for tempered stable 
and layered stable processes. 

In section 2 we begin by studying the Lamperti stable distributions, which are multivariate in- 
finitely divisible distributions with no Gaussian component and whose Levy measure is characterized 
by a triplet (a, /, ct), more precisely an index a S (0,2), a function /, and a finite measure cr, both 
defined on the unit sphere in . In particular, the radial component of any of these Levy measures is 
asymptotically equivalent to that of an stable distribution, with index a, near zero and has exponential 
decay at infinity. These distributions have a density with respect to the Lebesgue measure, have finite 
moments of all orders and exponential finite moments of some order. In the one dimensional case, the 
density is C°°. In some particular cases, we also prove that these distributions are self-decomposable. 

In section 3, we formally introduce the Lamperti stable processes and study their properties with 
emphasis in the one dimensional case, where we obtain an explicit closed form for the characteristic 
exponent. Motivated by the works of Rosihski [3T| and Houdre and Kawai [16], we prove in section 4, 
5 and 6 that Lamperti stable processes in a short time look like a stable process while in a large time 
scale they look like a Brownian motion, that they are absolute continuous with respect to its short 
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time limiting stable process and they admit a series representation that allows simulations of their 
paths, respectively. 

In the last section, we study some related processes: the Ornstein-Uhlenbeck processes whose 
limiting distribution is a Lamperti stable law and the Levy processes with no positive jumps whose 
descending ladder height process is a Lamperti stable subordinator. Finally we illustrate with several 
examples the presence of Lamperti stable distributions in recent literature. 



2 Lamperti stable distributions. 

In this section, we define Lamperti stable distributions on IR'^ and establish some of their basic 
properties. According to Theorem 14.3 in Sato |33j . the Levy measure H of a stable distribution with 
index a on IR"* in polar coordinates is of the form 

n(dr,dO = r--("+iMrcr(d^) 

where a € (0, 2) and ct is a finite measure on 5"*^^, the unit sphere on IR''. The measure a is uniquely 
determined by LI. Conversely, for any non-zero finite measure a on S'^^^ and for any a G (0,2) we 
can define an stable distribution with Levy measure defined as above. 

Motivated by the form of the Levy measure of the processes mentioned in the introduction and the 
previous discussion, we define a new family of infinitely divisible distributions that we call Lamperti 
stable. 

Definition 1. Let fj, be an infinitely divisible probability measure on JR"^ without Gaussian component. 
Then, fj, is called Lamperti stable if its Levy measure on IRq := IR'' \ {0} is given by 

C^(S)=/ a{dO rtBHV^^^He''-l)-^"+'^dr, B G 6(0, (2.1) 
Js''-^ Jo 

where a G (0,2), a is non-zero finite measure on S"^^^ , and f : S"^^^ ^ R is a measurable function 
such that 7 :— sup^^gd-i /(^) < a + 1. 

Note that v^'^ is indeed a Levy measure on IRq. To see this, we need to verify that 

(1 A \\x\\^)iy:^'^dx) = [ a{d^) H {1 A r2)e''-^(«)(e'" - l)-("+iMr < oo. 



On the one hand, since e''-^*^^-' (e'' — l)^("+i) ^ ,--("+1) as r ^ 0, [^we have that 

/ (T(dO / r2e'^-^(«)(e''-l)"^"+^Mr < 
JS'i-i Jo 

One the other hand, from elementary calculations we deduce 



oo. 



a(dO 

gd-l 



where 7 := sup^g^d-i /(^). Since 7 < a + 1, the above integral is finite and therefore is a Levy 
measure. 

In the one dimensional case / takes only two possible values, since 5*° = { — 1, 1}. In the sequel, 
we denote these two values by /(I) := /3 and /(—I) p ■ From the definition of the measure ct, we 
have ct({1}) = c+ and ct({— 1}) = c_. Therefore each distribution associated to the Levy processes 
mentioned in the introduction belongs to the class of Lamperti stable distribution. Following the 
notation of the introduction, for: 

• the subordinator 5, j3 — \ and c_ = 0, 

• the process /3 = 1 and p = a, 



^We say that / ~ g as a; ^ if for xq 6 M'' WvUx^xq f{x)/g{x) = 1- 
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• the process /? = ap + 1 and p = a{l — p), 

• the process (3 = ap and p = a{l — p) + 1, 



• the class of Levy processes with no positive jumps considered by Patie [2S], p = 1 — 1) and 
c+ = 0. 

Note that Laniperti stable distributions satisfy the divergence condition, i.e. 

/>oo 

/ e''^(^)(e''-l)-("+i)rfr = 00 for any ^ e S''-\ 
Jo 

Thus from Theorem 27.10 in |33], we deduce that they are absolutely continuous with respect to the 
Lebesgue measure. Also note that Lamperti stable distributions do not belong in general to the class 
of tempered a-stable distributions (see [S]), nor to the class of layered stable distributions (see [T^] 
for their definition). This follows respectively from the facts that the function 

e'-/(«) 

g(r £) = 7-^+" 

(e''-l)"+i 

is not completely monotone for any fixed ^ £ 5"^^^, and that for each ^ G S^^-^ 

IVtO -(a+l-f(n)r 

r-^ e as r — > cx). 

The above remarks motivate us to see Lamperti stable distributions as a particular class and study 
some of their properties. 

Proposition 1. Let p be a Lamperti stable distribution with Levy measure v^'^ given by (2.1), then 

WxW^ p{dx) < CO for all p > 0. 



Let 7 :— sup^ggd-i /{£,)■ ^/C < « + 1 ~ 7; then 



e'^l'^V(c^a;) < 00, 
In particular, for k < a + I and if f = k, we have 

e^"^"/i(c?a;) < 00 if and only if < a + 1 — k. 

Proof: According to Theorem 25.3 in Sato [33], we need to verify that the restriction of i^^'^ to the 
set {x € IRq : ||a;|| > 1} has the corresponding moment properties. To this end, we consider 

\\x\\Piy^'^idx) = f a{dO r rPe^(«''-(e'- - l)-("+i)dr 



which is finite since 7 < a + 1 . 

Now, we turn our attention to the exponential moments. Consider 



J{\\x\\>l} Js-i-^ Jl 

/>oo 



which is also finite for Q < a + 1 — 7. Next, we suppose that f = k. The former arguments imply that 
for ^ < a + 1 — K, the Lamperti stable distribution p has a finite exponential moment of order ^. In 
a similar way, it is clear that 

/•oo 

{||x||>l} Jl 
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This implies that / e'^^^''^^^i^^''^{dx) is finite if and only if ( < a + 1 — k. ■ 

J{\\x\\>l} 

Our next result shows that Lamperti stable distributions are self decomposable in some cases, and 
that they belong to the Jurek class, i.e. the class of infinitely divisible distributions for which the 
Levy measure v takes the following form 

/• /"OO 

KS)=/ <j{dO lBHm,r)dr, for B E B{Mi), 

Js-'-^ Jo 

where l{£,,r) is measurable ^ e S"^^^, and decreasing in r e (0,oo). 

Proposition 2. Let ji he an Lamperti- stable distribution on IR'' with Levy measure v^'-^ given by 



(2.1), then fi belong to the Jurek class. Moreover, /i is self decomposable if and only if f{^) < a + 1/2, 
for all^e S'^-^ and a £ (0,2). 

Proof: We first prove that Lamperti stable distributions belong to the Jurek class. To do so, it is 
enough to show that the function 

e/(«)'- 
^ (e'' ~ ' 

is measurable in ^ G S'^^^ and decreasing in r > 0. The measurability of ^ in ^ G S'^^^ is clear from 
its definition. To prove that £ is decreasing in r > 0, we fix ^ e S"^^^ and consider the derivative of 
£i(-)=£(^,-),i-e. 

Hence ([{r) < for r > 0, since /(^) < a + 1. This implies that ^i(r) is decreasing for r > and 
prove the statement. 

Now, we prove the second part of the Proposition. From Theorem 15.10 in Sato [33 we need to 
verify that the function 

'=(^''^)^(e^-l)"+i"' 

is measurable in ^ e S"^^^ and decreasing for r > 0, if and only if /(^) < a + 1/2 for all ^ G S"^^^ and 
a £ (0, 2). The measurability of fc in ^ G S"^^^ is clear from its definition. Fix ^ G S'^^^, and assume 
that /(f) < a + 1/2. Consider the derivative of h{r) = /c(f , r), 

h'{r) = _ ^^^^^ (e'- (1 + fiOr ^ {a + l)r) - /(Or - l) . 

Therefore, the function h is decreasing for r > if 

air) = e^il + (/(O -ia+ l))r) - /(e)r - 1, 
is negative for r > 0. Again, straightforward calculations give us that 

g'ir) - e'- (/(O ~a + /(e)r - (a + l)r) - /(f) and 

g"{r) = e'- (2/(0 - 2a - 1 + /(^r ~ (a + l)r) . 

It is not difficult to see that for /(O < a + 1/2, we have that g"{r) < for r > 0, which implies that 
g'{r) is decreasing for r > 0. Now, since ^'(0) = —a < 0, we deduce that the function g'{r) < for 
r > 0, which implies that g is decreasing for r > 0. Next, since ^(0) = 0, we have that the function g 
is negative for r > 0, which proves that h is decreasing for r > assuming that /(O < a + 1/2, for 
allfG S"'-!. 

Finally, it just remains to prove that if /(O > a + 1/2, for some f G S"^^^ and a G (0, 2), there exists 
a To > such that g^ro) > 0. This follows from similar arguments used above (take for instance 
/(O = a + & + 1/2, for 6 > 0, and a close to 0). ■ 

We finish this section with some properties of Lamperti stable distributions defined on IR. The 
first of which says in particular that the density of any Lamperti stable distribution belongs to C°° . 
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Proposition 3. Let fi be a Lamperti stable distribution on IR, then fi has a C°° density and all the 
derivatives of the density tend to as \x\ tends to oo. 

Proof: Recall that the function / takes two values, /? = /(I) and p = /(—I) as usual. According to 
Proposition 28.3 in Sato j33j and the form of the Levy measure of fi, it is enough to prove that 

g{r) = (-^^ip^d^' verifies that hm inf ^ > 0, (2.2) 

for some a G (0, 2). But this is immediate because for r sufficiently small, we have 



X 



p0x l-r 2 

2 dx> K -^dx = Kr^-", 



where K > Q. In particular when a = a, the condition in (2.2 1 is satisfied and the statement follows. I 



Before we state the last result of this section, we recall the definition of a particular class of 
distributions which is important in risk theory (see for instance fl2] and |21j'). 

Definition 2 (Class Take a parameter (5 > 0. We shall say that a distribution function G on 

[0,cx)) with tail G :~ 1 — G belongs to class L^^'^ if G{x) > for each x>0 and 

hm ^^^^ — — = p/^ for each x e R. 

u^oc G{u) 

The tail of any (Levy or other) measure, finite and non-zero on (a;o,oo) for some xq > 0, can be 
renormalised to be the tail of a distribution function and by extension, then is said to be in L'^^\ if 
the associated distribution function is in C'^^^ . 

Now, we will prove that the tail of the Levy measure of any Lamperti stable distributions defined 
in IR belongs to the class £("+1-/5)^ where (3 = /(I) as usual. 

Proposition 4. Let p, be a Lamperti- stable distribution on IR, then the tail of its Levy measure belongs 
to the class £("+1^/3)^ particular when /i is defined on IR_|_, we have that ji belongs to the class 

£(a+l-/3)^ 

Proof: First, we define 

roc ^f3r 

where K — I ^ — r^dr. Note that v corresponds to the distribution function associated to 

ii (e'--l)"+i 

the tail of the Levy measure of a Lamperti stable distribution. 
From elementary calculations, we get 

V{u-x)_r ( r dr' " 



v(u) = — / r r-rdr, M > 1, 



< (a + 1 - /3)e"("+i-« / . —^-rdr 

- V f^' I ^ (gr _ iy,+i 



-Q-1 



Similarly 



^^r^>(l-e-")"+i(a + l-/3)e"("+i-« e^P-(-+^)>dr 

I'iu) Ju-x 
— g(a+l-/3)a;^-[^ _ ^-u-^a+l ^ 

Therefore taking u large enough, we deduce that v e £("+1-/3). The case when /j, is defined in IR_|_ 
follows from Proposition 3.4 in Kyprianou et al. |21j . ■ 
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3 Lamperti stable Levy processes. 

Here, we introduce the class of Levy processes which is associated to Lamperti stable distributions. 
We also discuss, specially in the one dimensional case, a number of coarse and fine properties of their 
paths which are of particular interest for applications. 

Definition 3. A Levy process without gaussian component, and linear term 0, is called Lamperti 
stable with characteristics (a,f,a,9) if its Levy measure is given by (2.1). 

In the sequel, we denote the Lamperti stable Levy process with characteristics (a, /, cr, 6) by 
~ {X^,t > 0). Its characteristic exponent is defined by E[exp{i{y, Xj^))] = exp{—t'i'{y)) for 
t>0,y€M'^ where 

^{y)^z{y,e)+ [ (l-e'<^'->+z(2/,x)I{||,||<i})i.^-^(dx), (3.3) 



the measure v'^'^ has the form given in (2.1 1 and 6* G IR . 

The first property in study is the p-th variation of Lamperti stable processes. In particular, we 
prove that their p-th variation is similar to that of stable processes. 

Proposition 5. Let X^ be a Lamperti stable process with characteristics {a, f,a,d). 

i) Lf a E (1,2), the process X^ is a.s. of finite p-th variation in every finite interval if and only if 
pe (a,2). 

ii) The process X^ is a.s. of finite variation in every finite interval if and only if a E (0, 1). 

Proof: (i) From Theorem III in Bretagnolle [5], we have that for p E (1,2), the process X'" is a.s. of 
finite p-th variation on every finite interval if and only if 

\\x\\Pi'^'f{dx) < oo. 

{Ibll<l} 

Recall that 7 := sup^ggd-i /(^). From the form of the Levy measure i/"'^ and some elementary 
calculations, we have 



'{||x||<i} 
On the other hand, we have 



1 



\\xri.^^f{dx) < a{S''-')e'^ ___dr < a{S''-')e'^ r^-^'^+^Ur. (3.4) 



> e S^-^ : /(O > 0}) I ^^.^c 



-dr 



^ I{/(0<o}e^(«V(dx) (^^3I)^dr (3.5) 

> e S"-' : /(O > 0}) ^ I{/(«)<o}e^^«V(dx)) rP-("+i)dr, 

for some K > 0. Therefore X^ is of finite p-th variation on every finite interval if and only if p > a. 
The proof of part (ii) is very similar. According to Theorem 3 of Gikhman and Skorokhod |14j . it is 
enough to prove that 

/ \\x\w:'''{dx) < 00, 
-'{ii^ii<i} 

if and only if a £ (0, 1). But this follows from (3.4 1 and (3.51 taking p = 1, which concludes the proof. 



Recall that the characteristic exponent of a Levy process has a simpler expression when its sample 
paths have a.s. finite variation in every finite interval. In this case, we have that the characteristic 



exponent in (3.3 1 takes the form 



iy) = ~-z{d,y)+ f (l-e^<^'->)c^(da;), 
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where d ^ - / ..-/(d.) and, hr this case, d is known as the drift coefficient. 

-'{||:r||<i} 

In the rest of this section we work with real valued processes. We are now interested in two 
important properties of Levy processes: creeping and the regularity of 0. Define for each a; > 0, the 
first passage time 

T+ = inf |i > : > xj, 

with the convention inf = oo. We say that the Lamperti stable process creeps upwards if for all 
X >0, Po{X\ — x) > 0. If —X^ creeps upwards, we say that X^ creeps downwards. 

Proposition 6. Let X^ he a Lamperti stable process with characteristics {a, f,a,6). 

i) Lf a E (0, 1) and d > 0, the process X^ creeps upwards. 

ii) Lf a E (1,2) the process X^ does not creeps upwards. 

Proof: (i) The statement follows directly from part {i) of Theorem 8 in [25]. 
(m) From Theorem 8 part (Hi) in [22 , all we have to prove is: 

r xi^^-f{[x oo)) ^^^^^ ^^^^^^ H{x)= ( r iy^-f{{~oo,u])dudy. (3.6) 
Jo Hyx) J~xJ-i 



To this end, we study the behavior of the integrand in (3.6 1 near 0. We ffi'st note that 



ru -px 

/ 7 T^rrrrdr - 1 as u T 0, (3.7) 

which follows from the fact that e~^''(e^^ — l)"("+i) ~ x"^""''^-' when x increases to 0. Thus, it is 
not difficult to deduce that 

x"-^H{x) ^ ^ — as xlO. 

(2 — a)(a — l)a 



On the other hand, similar arguments as those used in (3.7) give us 



a 

Hence, 



"i/^'^([a;,c5o)) - - as xlO. (3.8) 



(2-a)(a-l)- as x i 0, 



which implies (3.6 1. The proof is now complete. 



Proposition 7. For a Lamperti stable process X^ with characteristics {a, f,a,9), the point is 
regular for (0,cxd) if one of these two conditions hold: 

i) ae [1,2). 

ii) a € (0, 1), and d > 0. 

Proof: Here, we apply Theorem 11 in [5?. When a e [1, 2), the process has unbounded variation 
(see Proposition 5 {ii)). Hence from part (i) of Theorem 11 in [21] we deduce that is regular for 
(0,oo). 

Now, we study the case where a £ (0, 1). In this case, the process X^ has bounded variation and 
according to part (ii) and (Hi) of Theorem 11 in [22J, the point is regular for (0, cx)) if the drift 
coefficient d > or if d = and the following condition holds 



^!|^M = oo, where H{x) - -y)dy. 
Jo H{x) Jo 



(3.9) 
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Thus, we only need to verify (3.9). In order to do so, we recall (3.8 1 

1 



2/"Cn(-oo,2/]) 



which implies that 



We observe then, that 



1 



a{l — a) 



U{x) 



which implies (3.9 1 and the proof is now complete. 



as y i 0, 
as a; J, 0. 



a(\ — a), as a; J, 0, 



Our next result deals with the computation of the characteristic exponents of Lamperti stable 
processes. Denote by 

r(z + a) 

which is known as the Pochhammer symbol. 

Theorem 1. Let he a Lamperti stable process with characteristics {a, f,a,9). 

i) If (X (0, 1) U (1,2), the characteristic exponent of is given by 

^lW = tXO - c+Ti-a) ii-iX + 1 (1 - /?)„) 

-c_T{-a){{t\ + l-p)^-{l-p)^), XeR. 

a) If a — 1, the characteristic exponent of X^ is given by 

^'i(A) = iX9 - c+ ({-iX + 1 - f3)iji-iX + 2 - /?) - (1 - (3)^(2 - (3) 



Where 6 is given by 



c^[{iX+l-p)i;{iX + 2-p)-{l-p)i;{2~p) , Ae 



-d if a e (0,1), 

- (c+a/3 - c_6p + (c+ - c_)(l - C)) ifa^l, 

c+ - c_ 



- c+ap - c^bp + 



a-1 



(3.10) 



z/aG (1,2), 



0,13, bp are given in (3.26), (3.21), respectively; C is the Euler constant, and is the Digamma function. 



Proof, i) First we will consider the case where a £ (0, 1). Without loss of generality, we assume that 
d = 0. Since a e (0, 1), we know that the characteristic exponent of is given by 



^l(A) = - c+/ {e'^^-l)-, -—^dx + c- (e*^"-l)^ T^r^dx 

^' \ .In ^ ^(e^-l)"+i .L^^ ^(e-^-l)"+i 



(e-^ - 1)^ 



(3.11) 



We compute each of these integrals which we call Ii and I2 respectively. Since all the computations 
involved are valid for all A G M, we center our attention in the variable /?. In order to compute Ii 
explicitly we will define in the set ?7 = {zgC:5R(z) <a+l}, the following function F : U C, 
given by 



F{z):^ (e^^"-l)- \, 

W \ ^(e--l)"+i 

»i 



e 



-zix 



dx^ I (e'^"-l),, 

^ ^(l-e-^)"+i 



dx 



(3.12) 
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where zi — a + I — z and ^{zi) > 0. Then by making an integration by parts in the last integral of 
(|3.12| we obtain for ^{zi) > 1 



a 



u'^-^{l - u)-°'du. 



(3.13) 



Now recalling the integral representation for the Beta function, (see we have for 5R(a),3?(5) > 

"'0 



we can express (3.131, in the following form: 



[u-'^ - l)u^i-i(l - = 







r(a + 6) 

-(a + zi - 1) T{-i\ + zi - l)r(l - a) 
a r(— iA + zi — a) 

(zi-i)r(zi-i)r(i-a) 



(3.14) 



a 



r(zi — a) 



finally by the recurrence relation for the Gamma function, r(a:; + 1) = xT{x), and the fact that 
zi = a + 1 — 0, we obtain 



F{z) = r(-«) 



r(-iA + a + l-z) r(a + l-z) 



r(-iA + 1 - z) 



r(i-z) 



(3.15) 



for 5R(zi) > 1, i.e. 5ft(z) < a. So we have the desired result for (3 < a. In order to obtain it for 
/3 e [a, a+1) we do the following: The equality (3.15) is valid in particular in Da = {z G C : ||z|| < a}, 
in order to extend it to the case where ||z|| < a + 1, we will prove first that both sides of the equality 
in (3.151 are analytic functions in the disk Da+i = {z G C : ||z|| < a + 1}. 
First we take F, and then using a series expansion we have 



F{z) 



(g^A^ _ 1) 



(e^ - l)"+i 



dx 



_ l)if^(e^ _ l)-("+i)da;. (3.16) 



n=0 



n 



Now consider the following 



E 

n=0 



< 



dx 



^ / 1 1 II \n 

y(|A|x)^^^(e--l)-("+i)dx 



n=0 



(|A|a;)ell"ll^(e^ - l)-("+i)da;, 



(3.17) 



which is finite when ||z|| < a + 1, therefore we can apply Fubini's Theorem in (|3.16|) and obtain 

OQ 



n=0 



(3.18) 



for z lE W, where 



-dx. 



which implies that F is analytic in Da+i- 

Since for ||z|| < a + 1 we have that 5R(— iA + a + 1 — z) > 0, and ^{a + 1 — z) > 0, therefore the 
function G : U ^ C, given by 



G'(z) =r(-a) 



r(-iA + a + l-z) r(a + l-z) 

r(-iA + 1 - z) r(i - z) 
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is analytic en Da+i- Since F and G are analytic in Dq+i, and _F = G in Da, we conclude that F — G 
in £>Q+i, which implies that 



h = F{/3) = T{-a) 



T{-iX + a+l- /3) T{a + 1-P) 
T{-iX + l-l3) T{l-(3) 



for all /3 < a + 1 . 

Now we compute the second integral in the right-hand side of (3.11 1 

1-0 



-px 



-dx 



(e-^ - l)"+i 

where pi = a + 1 — p, hence following the same arguments used in the computation of Ii , we get 

'^='^^~"^ I r(zA+i-p) — f(wr 

for all p < a + 1 . Therefore from the form of Ii and I2 we get 

^-lIA) - -c+T{-a) + 1 -/?)„- (1 - /3)„) 

- c_r(-a) ((jA + 1 - p)„ - (1 - . 

for all /?, jO < a + 1. 

Now we consider the case where a G (1,2). As in the case where a E (0, 1), we assume that 6 = 0. 
Since a £ (1, 2), the characteristic exponent of is given by 



*l(A) = -(^c+^ (e*^" - 1 - jAxI{,<i}) 



(e^ - l)«+i 



dx 



(e-^ - 1)" 



+1 



dx 



(3.19) 



We call Ii and /2 respectively the integrals in (3.191. To study Ji to do that we define the function 
G:U- 



given by 



G(z) := / (e*^- - 1 - zA.Tl{,<i}) 



(e^ - l)"+i 



dx 



(1 - e-^)"+i 



-da; -|- iA 



1 a;(e- 



+ iX e 



. 1 — x — e ^ 
(1 - e-^)"+i 



-f «A 



= lAa -|- iXI{z) + 
where zi = a + I — z, 5ft(zi) > 0, 



(1 - e-^)" 
^ {u-'^ - l)u^i-i -hiA(u- 1) 



(l~e-)"+i 
dx 



dx 



(l-u) 



a+l 



du. 



(3.20) 



e "^T- — ^—rdx- 



(1 - e-^) 



dx. 



1 (1-e-^)" ' 



and / : [/ ^ C is defined by 



J(z) := 



^ a;e-^(l - e-("-^)^) 
(1 - e-^)"+i 



da;. 



(3.21) 



We consider the last integral in (3.20), an integration by parts gives us for '^{zi) > 2 



^ (u^''^ - l)u^i-i -(-iA(u- 1) iA (;2i-l) 



{\-u) 



a+l 



du - 



a a 
1 



{u-''' - l)u^^-\l - u)-"du 



iX 



a 



(3.22) 
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We now compute the integrals in the right-hand side of (3.221, by making an integration by parts, 
rccaUing the integral form of the Beta function (3.14), and using the recursion formula for the Gamma 
function we obtain for the first integral the following 



1 - a 
_ {^iX + zi - 2) 
I — a 

= T{l-a) 



iX 

1 - a 







^ 



1 - a 

rHA + zi) r(zi) 

(— lA + zi — l)r(— iA + zi — a) (zi — l)r(zi — a) 



(3.23) 



and for the second integral in the right-hand side of (3.221 
»i 

{u-"^+''-^ -l){l-u)-°'du 



J 

Jo 



1 (-iA + zi-2) 



-iX-\-zi—3 



(1 



1 (-iA + zi - 2) r(-a + zi - 2)r(2 - a) 

a — I 1 — a r(— iA + zi — a) 

1 r(-iA + zi) r(l-a) 

a — I {—iX + zi — 1) r(— iA + zi — a) ' 



(3.24) 



so using (3.231 and (3.24), in (3.22) and recalling that zi = a -I- 1 — z, we get 



^ (u-'^ - l)u^^-^ + iX(u - 1) , iX , 
(1 ~ a — 1 



r(-iA + a+ 1 - z) r(a+l-z) 

r(-iA + 1 - z) r(i - z) 



if 5R(z) < a - 1. 

So if we consider the function P : U 



P{z) 



iX 

a-l 



+ r(-a) 



defined by 

r(-iA + a + l-z) r(a + l-z) 



(3.25) 



r(-a + i-z) r(i-z) 

we have that 

G'(z) = iAa-HiA/(z)-|-P(z), 

for 3fi(z) < a — 1, in particular the equality holds in the set -Dq_i = {z G C : ||z|| < a — 1}. 

We will prove that G is analytic in -Da+i, so we follow the same method as in (3.17) and (3.18) and 

obtain 

°° poo 

Gi^) = Y ^ (e'^" - 1 - iXxl!,<n) a;"(e" - iy^"+^^dx 

oo 

= ^a„z", 



n=0 



for z e I?Q+i, where 



In a similar way we prove that the function / defined in ( |3.21 1 is an entire function. 

We note that if z e -Dq+i, then 3fJ(— lA -I- a -I- 1 — z) > 0, and 5R(q; -|- 1 — z) > 0, which implies that the 

function P defined in (3.25) is also analytic in D^+i- 

Finally since G — iXa + iXI + P in Da-i, and both sides are analytic in Da+i, then G — iXa + iXI + P 
in Da+i. 

Therefore the first integral in the right-hand side of (3.19), Ii, is given by 

1 



h = G{(3) = iX[a[3 + 



a-l 



+ T{-a) 



T{-iX + a + l- (3) T{a+1-I3) 
T{-iX + 1-/3) r(l-/3) 
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where 

— a + I{(3) 



(1 - e-^)"+i 



dx 



- -da; 



-dx, (3.26) 



for all /3 < a + 1. Now we compute the second integral m the right-hand side of (3.19 1 

{u'^ - l)uP^-^ + i\{u - 1) 



-i\b„ 



(1 



-du, 



where 



(1 - e-^)"+i 



da; - 



1 — a; — e 



-da; 



(1 - 6"=^)" 



-da;, 



(l-e-^)°+i 

and pi = a -|- 1 — p, hence following the same arguments used in the computation of Ii , we get 

1 



(3.27) 



h = -iX bp + 



a - 1 



T{i\ + a+l- p) T{a + l-p) 

r{tx + i-p) r(i - p) 



for all p < a -|- 1 . Therefore from the form of Ii and I2 we get 



- c+T{-a) ii-iX + 1 -/?)„- (1 - /3)„) 



a-1 

- c_r(-a) {{iX + 1 - p), - (1 - p)„) , 
for all (3, p < a + 1. 

ii) Now we will compute the characteristic exponent when a = 1. In the following we assume that 
9 = and that c+ = c_ = 1. Since a = 1, the characteristic exponent of is given by 



*l(A) 



{a=<l}J 



(e^ - 1)' 



rda; 



iXx 



- 1 - «AxI{:r>_i}) 



rda; 



(e-^ - 1)2 

We will follow the same arguments used in the first part of the computation of the characteristic 
exponent in the case a G (1, 2). But to compute the two integrals in (3.22) we will need the following 
integral representation for the Digamma function (see |15j ) 



' V'-^ - 1 
z - 1 



dt - C, for z e C, 



(3.28) 



where C is the Euler constant. 

Now by making an integration by parts, using (3.281, and the recurrence relation for the Digamma 
function ■(/'(z -|- 1) = i'iz) + z^^, we can express the first integral, for 5R(zi) > 1, in the following form 



f\u-'^ ~ l)u'^~^{l - u)-^du = - (^J^ 



-du 



u-l 

= V(2i " 1) - V^-iA + Z1-I) 
1 



1 _ 1 



U- 1 



-du 



zi-l 



ipi-iX + zi) + 



~iX + zi — 1' 



(3.29) 



As for the second integral in the right-hand side of (3.22) 



(y-»A+2i-2 _ ^^1^^ _ u)-^du = -i;{~iX zi - 1) - C 







1 

-lA -I- zi — 1 



— 'ip{~iX + zi) — C. 



(3.30) 
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So using (3.291 and (3.301, in (3.221 and recalling that zi = 2 — z, we get 



(1-U)2 



du = iX{l -C) + {~iX + 1 - z)i/;(-iA + 2 - z) ~ (1 - z)V'(2 - z) 

(3.31) 



if 5R(z) < 1. 

We note that (3.31 ) can be extended to the case where, 5R(z) < 2, by the same arguments used in the 
case a G (1,2), we only need to remark that the function P : U ^ C defined by 

P{z) ^ iX{l -C) + {-iX + 1 - z)^j{-iX + 2 - z) - (1 - /3)?A(2 - z), 



is analytic in the disc D2 = {z e C : ||z|| < 2}. This implies that (3.311 is true for all z G D2, so in 
particular 



Ji = G{(3) = iX{dp + 1 - C) + {-iX + 1 - /3)^(-zA + 2 - /3) - (1 - /3)V'(2 - /3), 



where 



1 a;e-'^(l -e-(i-^)^) 
(1-6-^)2 



I "e-V^"^' A (1-e-) 



for all /3 < 2. 

Finally by the same arguments used in the computation of Ii, we obtain 



/2 = -iA(6p + 1 - C) + (iA + 1 - p)^{iX + 2 - p) - (1 - p)V'(2 - p), 



where 



^ a;e-^(l - e^^i-'')-^) 



/o (l-e-)2 ,0 
for all p < 2. Therefore from the form of Ii and I2 we get 



/ (i-e-^)2 '^'' + / 



«'l(A) = -iA [^c+dp - c^bp + (c+ - c_)(l - C) 

- c+ ( (-zA + 1 - PM-tX + 2-P)-{l- /3)^{2 - 13) 



c_ (iA + 1 - p)ij{iX + 2 - p) - (1 - p)^j{2 - p) 



for all /3, p < 2. 



Using the well known relationship between the Laplace and the characteristic exponents, we obtain: 

Corollary 1. Let he a Lamperti stable process with characteristics {a, f,<J,d). 

i) Let a G (0, 1) and suppose that is a Lamperti stable subordinator, then its Laplace exponent 
is given by 

$l(A) - dA - c+T{-a) f (A + 1 - /3), - (1 - (3)^ , A > 0, 



where d > 0. 

ii) Let a £ (1,2) and suppose that X^ has no positive jumps then its Laplace exponent is given by 
$i(z) = -^A + c_r(-a) f(A + 1 - p)„ - (1 - p)^) , A > 0, 



where 6 is given by {3.10). 

Remark 1. This Corollary has, as particular cases, the two recent results found in: 

i) Corollary 2, and Lemma 4-, UOl, where the result is obtained by means of the Lamperti 
transformation. 
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a) Proposition 3.1 in J28l, where the Laplace exponent is obtained using special functions. The ideas 
in [SBJ as well as in 161, inspired parts of the proof if Theorem 1, specifically the decomposition 



(3.20) 



Now we turn our attention to other group of properties. Let H = iHt,t > 0) be the increasing 
ladder height process of X'" (see chapter VI in [2 ) and H = {Ht,t > 0), its decreasing ladder height 
process. Denote by k and k for the characteristic exponents of H and H, which are subordinators, 
and suppose that drifts to — oo and i^"'-''(0, oo) > 0. Under this hypothesis, the process if is a 
killed subordinator and we denote by Hh for its Levy measure. The following result give us a relation 
between i^^'^ and Hh. 

Proposition 8. Let X^ be a Lamperti stable process with positive jumps and characteristics (a, /) 
such that it drifts to — oo. Then, the tail of the Levy measure of H , its increasing ladder height process, 
belongs to £("+i~^) and 

ly^'-^ (u, oo) ^ k{—i {a + 1 ^ (3))Il{u,oo) as u — > oo. 

Proof: The proof follows directly from Proposition 5.3 in '21] and proposition 4. ■ 

We finish this section with some properties of Lamperti stable processes with no positive jumps. 
Proposition 9. Let X^ be a Lamperti stable process with no positive jumps and characteristics 



(a, p, a, 6), such that 9 ^ in (3.10). Then 



i) there exist po G (1,2) such that Xls drifts to oo, oscillates or drifts to —oo according as p ^ 
(-oo,po), Pa or p£ {po,a + l). 

ii) for p e (po, OL+ 1), we have that there exist A > such that 

Vo{s^^ > a;) - as x oo, (3.32) 

where = supj>qX/', c = -logPo(i?i < oo), k = Eo(i?ie^^i ; iJi < oo) and H is the 
increasing ladder height process. 

Hi) for p e (po, 0L+ \), we have that there exist A > such that 

Po(^I{X^) > x^ Kx-^, as x->oo, (3.33) 
where K is a positive constant and 

/(X^) = ^ exp{x,^}dt. 

iv) the process X^ has increase times. 

v) the process X^ satisfies the Spitzer's condition at oo, i.e. 



lim 

t — *oo 



1 r* 

- P(Xf > 0)ds = l/a as x oo. 
i Jo 

vi) the process X^ satisfies the following law of the iterated logarithm 



r X,^<i>l\t-Hog\logt\) _ 

nmsup — j — a [l — a) a.s., (3.34) 



x^o log I log 

where <i>^^ denotes the right-continuous inverse o/<i>^^ 
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Proof: (i) We know that in this case a E (1,2), so from Corollary VII. 2 in [5], the process drifts 
to H-oo, oscillates or drifts to — oo according as <i>'^(0'^) is positive, zero or negative. Hence, from the 
Laplace exponent of we have, using the recursion formula for the Gamma and Digamma functions, 
the following 

a>i^(0+) - c_r(-a)(l + a- p)„(^A(l - p+ a) - V(l - p)), 
T{l + a-p) 



c_r(-a)- 
-9ip)- 



r(3 - p) 



(2 - p)(l - p)ai,{l -p + a)-i,{l-p))+3-2p 



(3.35) 



We have from (3.35) that g{l) < 0, and g(2) > 0. On the other hand, in the interval (1,2), the 
function g is continuous and decreasing which implies that there exist po € (1, 2) such that g{po) = 0. 
Thus, we deduce that X^ drifts to oo, oscillates or drifts to — oo according as p e (— oo,po), P — Po 
or p g (po,a+ 1). 

(ii) First, note that any Levy process with no positive jumps which drifts to — oo satisfies that its 
Laplace exponent has a strictly positive root. Hence for a Lamperti stable process with no positive 
jumps and with p G {po, a + 1), there exists A > such that 

Eo(exp{AXi^}) =1, 

i.e. that X^ satisfies the Cramer condition. Thus, the main result in [3^ give us the sharp estimate in 
([332]). 

(Hi) First note that X^ is not arithmetic and that under our assumptions the Cramer condition is 
satisfied for some A > 0. Hence from Lemma 4 in f5^, we get the sharp estimate (3.33) for the 
exponential functional /(X^). 

(iv) Here, we need the following estimate of the Pochhammer symbol (see for instance [26J), 



(A+1 — p)q'^A" as A ^ oo. 
From Corollary VII. 9 and Proposition VII. 10 in [2] we know that X'" has increase times if 

/oo 
A-3$^(A)dA < oo, 



(3.36) 



which in our case is satisfied since from (3.36), we have 

$l(A) ~ c_r(-a)A" 



as A 



oo. 



(3.37) 



(v) From (3.37), we see that <I>l is regularly varying at oo with index a. Hence, the statement follows 
from Proposition VII. 6 in [^. 

(vi) Since is regularly varying at oo with index a, we have that its right-continuous inverse <i>^^ is 
regularly varying at oo with index 1 /a which corresponds to the Laplace exponent of the first passage 
time of X-^ (which is a subordinator). Therefore, from Theorem HI. 11 in [2] we deduce the law of the 
iterated logarithm (3.34). ■ 



4 Short and long time behaviour. 

Motivated by the works of Rosihski [31j and Houdre and Kawai [TB] , we study the short and long time 
behavior of Lamperti stable processes. In particular, we will show that this class of processes share 
with the tempered and layered stable processes, the peculiarity that in short time they behave like 
stable processes. 

The convergence in distribution of processes, considered in this section, is in the functional sense, 
i.e. in the sense of the weak convergence of the laws of the processes on the Skorokhod space and will 

be denoted by " " . 
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Proposition 10. Let be a Lamperti stable process with characteristics (a,/, cr, 0) and 



5d-i 
Then, 

(h-'^^" {Xf^^ - htri^),t > Oj ^ {Xt,t>0) as h^O, 
where {Xt,t > 0) is a stable process of index a. 

Proof: The proof is similar to that of the short time behaviour of layered stable process, since for each 

e/(?)'-(g'- _ _ ^-("+1) as r ^ 0. 

Thus, we follow the proof of Theorem 3.1 in [1^ with 

g(^,r) = e-''(«)''(e'^-l)^("+i) 
and the desired result is obtained. ■ 

Theorem 2. Let X^" be a Lamperti stable process with characteristics (a, /, cr, 0) and 



Va = 

Then, 



gd-i 



(^h-'^/^{Xl^t-htr]a),t>0^ S {Wt,t>0) as h ^ oo, (4.38) 
where (Wt,t > 0) is a centered Brownian motion with covariance matrix 

Proof: According to a standard result on the convergence of processes with independent increments 
due to Skorokhod (see for instance Theorem 15.17 of Kallenberg [H]), the functional convergence 



(|438| holds if and only if 

h-^^^ {X^ - hria) ^Wi as h ^ oo. 
Now, we introduce the following transform for positive measures, for any r > 

{Tri^){B) = iy{r-^B) for S G S(]R'*). 

Note that the random variable h^^/'^Xf^^ is infinitely divisible and since it has finite first moment, 
we may rewrite its characteristic exponent as follows; 

ih [ (y,x)I{|i,|i>i}(T,-i/.<'0(dx) - /^ / (e'<^^^> - 1 - x)I{|i,,i<i}) (T;,-V2i^^'0(dx). 

Hence, from Theorem 15.14 of Kallenberg [T^ we only need to check the following convergences as h 
increases: 

a) h(Tf^-i/2i'^'^) converges vaguely towards on IRq, 

b) for each fc > 0, hj a;a;'(r^-i/2i/"'^)(dx) ^ / xxV^'^(dx), 

J||2;||<K JR^ 



c) for each fc > 0, h x{T,,-i/2iy^'f){dx) 0. 

J\\x\\>k 
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We first prove (a) or equivalcntly 



lim / g{x)h{T,,-^,.v':-f){dx) = Q 



(4.39) 



for all bounded continuous functions g : IRq IR vanishing in a neighborhood of the origin. Let g 
be such a function satisfying that \g\ < C, and that for some 6 > 0, g{x) = on {a; G IRg : \\x\\ < S}. 
Let 7 := sup^^gd-i f{S,), then we have 



h gix){T^-.,.iy:-f)idx) 



<h'+^/' [ <j{dO r\gH)\e 
Js''-^ Jo 

aidoTlgirOl 3 



^^^(e'-''^^^-l)-("+i)dr 



Sd-i 



(4.40) 



On the other hand, since 7 < a + 1 it follows 

lim r^- — — = 0, 

rwoo [e^ — 1)"+1 

then for e > sufficiently small, there exist M > such that for all r > M 



Since r > 5, we may take h > {^)'^ in (4.401 and obtain 



g{x) {T^-^^.,y:'f){dx) 



IRS 



<e[ a{dO r \g{rO\^dr 



< eC 



Sd-i 



Note that the last integral in the right-hand side of the above inequality is finite and therefore the 
convergence (4.391 follows. 

Next, we prove part (6). First note that / ||a::|pt'^^''^(da;) is finite. This follows by similar arguments 

as those used in proposition 1. This implies that the integral / xx'v^'-^ (dx) is well defined. Now 
take fc > fixed, and note that 

h [ xx'{Tf,-i,2V^'f){dx) = [ xx'v^'f{dx) [ xx'iy^'f{dx), 

■^{\\x\\<k} J{\\x\\<h^/^k} JRg 

as h goes to 00, which proves part (6). 

Finally, we consider fc > and recall that 7 = sup^^gd-i f{£,), then 



{||x||>fe} 



< (1 )-("+!) 



^^^(e'-''^^^-l)-("+i)dr 



h^+^/^ / ^a{dO I re'-''''' f-^- dr 



-kh^'^{a+l-l) 



hk 



(1 _ f,-kh^/^^a+i \a+l-j (a + 1 - 7)^ 
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which goes to as ft, ^ oo since 7 < a + 1 . This completes the proof. 



Let us apply the above results to the special cases treated in the introduction. In particular 
when we start with a stable process {X,Fx),x > 0, of index a, applying the result in short time 
behavior after various transformations we return to this initial process. Recall that associated to the 
stable process three Levy processes are obtained via the Lamperti representation of pssMp: ^* ,£,\£,^- 
Then the normalization of any of them according to proposition 10, converges weakly in the space of 
Skorokhod to the original stable process X, i.e. 

^ km LT^ norrn ^ ^ 

X km X* ^ X'' ^ X"^ "'^^ Xt ^ X as h^Q 

where kill, LT, DT and norm means killing , the Lamperti representation of pssMp, Doob-transform 
or conditioning, and normalization of a given process, respectively. Moreover X'~^ is the conditioned 
process (to be positive or to hit continuously), stands for any of the Lamperti stable processes 
and £,* , and Xj^ is the normalization of each of them given in proposition 10. In the same spirit 
we could also write, using theorem 2, 

kill -t^^ LT norm d , 

X — > X — > X — > X^ — > as ft ^ 00, 

^ kill „^ DT LT norm „^ d , 

X — > X — > X — > X — > X,^ W as h ^ 00, 
where is a centered brownian motion. 

The final result of this section follows the line of reasoning of last remark but uses additional tools 
that we shall briefly introduce. In [7 the convergence in the Skorokhod space is studied in relation 
to the second Lamperti transformation (LT2), i.e. the one that transforms Levy processes with no 
negative jumps to continuous state branching processes. The problem of explosions is difficult to 
handle in this metric so the authors consider another metric doo on the Skorohod space which given 

by 

rfoo(/,5)-lA inf ||/-.goA||ooV||A-/||oo. 

and where Aoo is the set of increasing homeomorphisms of [0, 00) into itself. According to the authors 
the convergence in this metric implies it in the usual Skorohod metric. Two of the main results in 
[7] say that the Lamperti transform LT2, is continuous with this new metric (proposition 4) and that 
a sequence (K*'") of stopped Levy processes with no negative jumps converges in this new metric 
towards Y, a stopped Levy process when the sequence of the associated Laplace exponents of (y*'") 
converges towards the associated Laplace exponent of Y (proposition 5). Therefore, a combination of 
the results mentioned above and proposition 8 give us the following corollary. 

Corollary 2. Let X and X^ be a stable proces of index a with no negative jumps and a Lamperti 
stable processes with no negative jumps with characteristics (a,/?) which does not drift towards +00, 
respectively. Let Yh = LT2{Xl^) and Y = LT2{X). Then 

Yh^Y as h^O. 

5 Absolute continuity with respect to stable processes 

We showed that in short time a Lamperti-stable process behaves like a stable process, now following 
Rosihski [Slj we will relate the law of both processes. In other words, we will find a probability 
measure under which the law of a Lamperti stable process with characteristics (a, /, a) is the same 
that the law of the short time limiting stable process with index a. 

Theorem 3. Let P and Q be two probability measures on (fl,J-') and such that under P the canonical 
process {Xt,t > 0) is a Lamperti stable process with characteristics {a,f,a,a), while under Q it is a 
stable process with index a with linear term b. Let {T-t) be the canonical filtration, and assume that 
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i) P|jFj and Q\j^t are mutually absolutely continuous for every t > if and only if 



Sd-i 



if a & (0,1), 



a - < 



S<i-1 Jo 



Sd-1 



-(a+i)d. 



i/ae (1,2). 



For each t > 0, 



dQ 
dP 



JJt 



where (Ut,t > 0) is a Levy process defined on (fl,J-,P) by 



Ut — lini 

40 



E 



{5-e(0,t]:||AX,||>e} 



||AX,||/(AX,)(-p||AX 



l)-("+i)||AX,||"+i 



i(C^-n)({zeMM|z|| >e}) 



/n f/ie above right hand side, the convergence holds P-a.s. uniformly in t on every interval of 
positive length. 

Proof: From Theorem 33.2 in Sato [33], we only need to verify that 

l^^<p{x)/2 „ l)2n(da;) < oo. 



where (p : IRq ^ IR is defined by 



dz^ 

dn 



In particular, we have ^{rS,) = log [e^^'^^^e'^ — 1) ("+i)7'"+i^. Thus, we need to check 



5c!-l 



a(dC) 



n 2 



(e'- - l)("+i) 



- 1 



1 



dr < 



(5.41) 



By Taylor expansion and the Lagrange form for the residual, we have (e*" — 1) = re'"^'', where Or G (0, 1). 
This implies 

r -^e'-W^^-^-^^+i)). (5.42) 

Now, noting that /(O - (a + 1) < /(O - Ov{a + 1) < /(O, it follows 

g'-(/(0-(" + l))/2 _ I < gK/(«)-e.(a+l))/2 _ ;^ < gr/(C)/2 _ ^ 



and since /(^) < 7 = sup^gg.d-i /($), we have 



(5.43) 



Using a Taylor expansion again and (5.43), it is clear that there exists a constant i? > such that if 
r < i?, then 

„2 



(e'-(/(«)-«'-("+i))/2 _ 1) < i^(/2(0 + l)r2 



(5.44) 
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where K is a, positive constant. Hence from (5.42) and (5.441, it follows that 

2 



a(dO 



R 



1/2 



<K(a{S''-^)+ [ 



1 



-dr 



R J.2 



„l+a 



-dr, 



which is finite because a € (0, 2) and / e L (S 



2/Qd-l n/Qd-l 



), ct). In the case when r > R, we have 



Sd-i 



a(dO 



e'-/(«)r(i+") 



1/2 



n 2 



- 1 



-dr 



< 4 ( (1 - e"«)-("+i) / (T(de) / e'^(/«)-("+i»dr + a(5'^-i) 

/"OO 



-dr 



< MS""-') (1 - e-^)-("+i) / e'-('^-("+i»dr + 



R 



-dr 



which is also finite because 7 < a -|- 1. Therefore (5.41 1 follows. 



The proof of the second statement of the Theorem follows directly from Theorem 33.2 of Sato 



Here, we follow the same notation as in Theorem 3. Note that under the conditions of Theorem 4.1 
in [TB] , if i? is another probability measure on (O,^) under which the canonical process X = {Xt,t > 0) 
is a layered stable process, we have that R\jr^ and QIjFj are mutually absolutely continuous for every 
t > 0. From our previous result, we obtain the corresponding result for Lamperti stable processes, i.e. 
that -RIjTj and are mutually absolutely continuous for every t > 0. Similar result holds for the 
tempered stable processes, see Theorem 4.1 in [3T| . 

6 Series representations of Lamperti stable process 

In this section, we establish a series representation for Lamperti stable processes which allow us to 
generate some of their sample paths. To this end, we will use the LePage's method found in |22. We 
first introduce the following sequences of mutually independent random variables defined in [0, T] . Let 
{ri}i>i be a sequence of of partial sums of iid standard exponential random variables, {Ui}i>i be a 
sequence of uniform random variables on [0, T], and let {Vi}i>i be a sequence of iid random variables 
in 5''^^ with common distribution a{d£,)/a{S'^^^). In order to use the LePage's method, we consider 
the following function : (0, 00) x S"*"^ IR+ given by 

P^^iu,C) := inf > : p{[x,oo),^) < uj, 

where 

/■oo 

p([x,oo),0= / e-^«)'^(e'- - l)-("+i)dr. 
Now, let {ci}i>i be a sequence of constants defined as follows, 

C; = ^ ^E(p-l(s/T, Vi)ViI{^-i(,/T,y,)<i})ds. 

Then from Theorem 5.1 in [35], the process 

00 \ 

J2 {p-\T^/T,v.)v.l^u^<t} -c..-),te [o,T] , 
i=l ) 

converges uniformly a.s. towards a Lamperti stable process with characteristics (a, /, cr) and linear 
term ^? = (in the Levy-Khintchine formula). In particular when /(^) = 1, we have that 

p-i(u,e)=ln(l + (a?.)-i/"), 



21 



hence the series representation for a Lamperti stable Levy process with characteristics {a,l), is 
as follows 



where 



e(^i) / ^ + ) I{ln(l + (a.T-i)-i/=)<l}ds. 



Let us observe below some sample paths of this particular Lamperti stable process generated via the 
series representation. 



7 Associated processes and examples. 

Here, we are interested in study some related processes to Lamperti stable distributions (or processes) 
and give some examples of Lamperti stable processes which appear in the literature but they are 
not the main objects in study. In particular, we study the Ornstcin-Uhlenbeck process and the 
self-similar additive process related to a Lamperti stable distribution in the case when the latter is 
self-decomposable. We also investigate the parent process of a Lamperti stable subordinator. 



7.1 Ornstein-Uhlenbeck type processes and self-simlar additive processes. 

Ornstein-Uhlenbeck type processes appear in many areas of science, for instance in physics, biology 
and mathematical finance. One of the particularity of these processes is that its limiting distribution 
is self-decomposable. Recall that a random variable Y on IR'', distributed as a Lamperti stable law 
with characteristics {a,f,a) is self-decomposable if and only if / < a + 1/2. Therefore, according 
to Wolfe [35] and Jurek and Vervaat [18], there exists a Levy process Z = {Zt,t > 0) on IR'', with 
Eo(log+ |Zi|) < oo such that 

/"OO 







where c > 0. Consequently, one can define an Ornstein-Uhlenbeck type process driven by Z with 
initial state Uq and parameter c > 0, that is the solution of 

Ut = Uo + Zt-c f U,ds, 
Jo 

and such that the law of Ut converge towards the law of 1" as i goes to oo. From Theorem 17.5 in 
|33j . we have that the process Z has no Gaussian component, its Levy measure is given by 

/• /"OO 

J5<i-i Jo 

where 

MO = (g/!'|)l+. {re^im) ~a-l) + e^- rfiO l) ; 

and linear term 

V = cq — I - — -Iiz{dx). 



In the one dimensional case, the form of the Levy measure of Z is reduced to 

T\z{dx) = c(^c+^-^-^^^^p^ (a;/3 1 - + xe"(a -f 1 - /?)) I{.>o} 



xp—l + e ^ + xe "^(a + 1 — /?) j I{a;<o} Idx, 



where /3 = /(I), p = /(—I), c+ = cr({l}) and c_ = a{{—l}), as usual. In this case, we have another 
process which is related to the Lamperti stable distribution Y, to Z and to the Ornstein-Uhlenbeck 
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Figure 1: a = 0.5, / = 1, = a{-l) = 



1. 



v'f 



f 



1 2 3 4 5 



Figure 2: a = 1.5, / = 1, = a{-l) 



1 2 3 4 5 



Figure 3: a = 1, / = 1, a{l) = cr(-l) = ] 
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Figure 4: a = 0.5, / = 1, a{l) = 1, £7(-l) 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 5: a = 1.5, / = 1, C7(l) = 0, (t(-1) 



8 



5 10 15 20 



Figure 6: a = 1.9, / = 1, C7(l) = 1, a{-l) 
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type process U . To this end, recall that in Theorem 16.1 of [33]) it is showed that a distribution is 
self-decomposable if and only if for any fixed > 0, it is the distribution of Vi for some additive 
process V = {Vt^t > 0) which is self-similar, meaning that for each fc > 



We remark that self-similar additive processes can be used to model space-time scaling random phe- 
nomena that can be observed in many areas of science. In particular, they are recently used to model 
asset prices and the risk-neutral process (see for instance [5]) in financial mathematics. 

Assume that V is the self-similar additive process associated to /, in which case Vi has the same 
law /, and that H — c. From Theorem 1 in |17] . there are two independent copies of Z denoted by 
Z(-) = {Z'f \t> 0) and = {z'l^\t> 0) which are defined by 

(_) m (+) (ao. r' ^ 



The process V can be recovered by 



Vr 



CO 

e-^hizl'^ if < r < 1, 

log(l/r) 

/•log(i-) 

y + / e"*dZ}+' if r > 1, 
Jo 



and moreover = e-^^Vet.t > 0) is the Ornstein-Uhlenbeck process driven by Z^"*"^ with initial 

state / and parameter c; and = e*'^Vf,-t,t > 0) is the Ornstein-Uhlenbeck process driven by 

— with initial state / and parameter — c. 



7.2 Parent process 

Motivated in generating new examples of scale functions, Kyprianou and Rivero [23] constructed 
Levy processes with no positive jumps around a given possibly killed subordinator which plays the 
role of the descending ladder height process. One of our aims is to determine the characteristics of 
the Levy process with no positive jumps whose descending ladder height process is a Lamperti stable 
subordinator. 

Let be a Lamperti stable subordinator with characteristics (a, /3, cr, 9) with zero drift and no 
killing rate. Since the density of its Levy measure is decreasing, then according to Theorem 1 in |23j , 
there is = {Xf^^t > 0), a Levy process with no positive jumps that we call the parent process 

of X^ whose Laplace exponent is given by 

^pl{\) = X^l{X). for A>0, 

where is the Laplace exponent of X^. Moreover, the process X^^ has no Gaussian coefficient, its 
Levy measure is given by 

UpL{dx) = c+j-^^—^^{^{a + l- I3)e-'' + (3^dx for x <Q, 

and with linear term 

b= [ xIlpLidx). 

J(-oo,l) 

Note that oscillates or drifts to oo according to whether is equal zero or strictly positive. 

From the form of its Levy measure, we deduce that X^^ is the sum of two Lamperti stable processes 
with no positive jumps X^ and X^ with characteristics (a + l,/3 + l,cri,6i) and (a -I- 1, /?, cr2, ^2), 
where ai({l}) = (72({1}) = 0, 

ai({-l}) ^ c+{a + 1-/3), ^2({-l}) = C+/3, 

bi= .Till (da;) - a/3+1, h ^ xll2{dx) - a/s, 

J(-oo,-l) J(-oo,-l) 
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and Hi and 112 are the respective Levy measures of and X^. On the other hand, the binomial 
expansion give us 



J3x 



^0 n!(a + 1~/3)„ ' ^''^^^ 

which is clearly log-convex on (0, oo) since it is completely monotone. Hence according to Theorem 
in 2 j23], there is a subordinator X*'^ with Laplace exponent such that 

W)-^^ for A>0. 

Moreover the subordinator X*''" has no drift and no killing term and the scale function of the parent 
process X^^ is determined by 



W{x)^ / Ul{y,oo)dy, 
Jo 

where is the Levy measure of X*'^. Note that for P — 1, we have that 

Ill{y,^) = —-^ r(l-e-^r-\ for y > 0, 

1 (a)l (1 — a) 

but for /3 7^ 1 we do not have an explicit form for . 

The Example 2 in [23] is related to the Lamperti stable subordinators considered above but with 
a given killing rate. Let us explain in detail such example in terms of Lamperti stable processes. Take 
X^'^ to be a Lamperti stable subordinator with characteristics (a, /3, cr, 9) with zero drift and killing 
rate given by 

r(i-/3) 

According to Kyprianou and Rivero |23j . there is a subordinator, here now denoted by Y with no 
drift, no killing rate and Levy measure given by 

which is the sum of two subordinators, one of which is a Lamperti stable with characteristics (1 — 
a, /? + 1 — a, (Ty), where 

y and ay{{-l}) = Q. 

Moreover, the Laplace exponent of the subordinator Y satisfies that 

0f(A)=-4yy, for A>0. 

From the form of Hy, we have the restriction that (3 < \. Thus, his parent process Y^ , a spectrally 
negative Levy process, has Laplace exponent 

A^r(l-/3 + A) 
^^^(^)^r(l™/3+A + a)' 

which has no Gaussian component and its Levy measure satisfies 

Hyp (-oo,y) = T\Y{dy)/dy. 



According to Kyprianou and Rivero and by (7.451, its associated scale function is given by 



^ n!(a + 2-/3)„ V 7' " 



ri=0 
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Now Y*'^, the parent process of the Lamperti subordinator X^'^ with killing rate K, is a spectrally 
Levy process which drifts to oo, with Laplace exponent 

c+r(-a)Ar(A + l-/3 + a) 

^ nx+i-p) ' 

which has no Gaussian coefficient and whose Levy measure satisfies 

Hyp,, (da;) = IfpL(da;), x < 0, 

with linear term 

xIlpL{dx) - K, 

and the associated scale function is given by 

c+r^(f - a) Jo 

As the process X^^ , Y*'^ may be seen as the sum of two Lamperti stable process with no positive 
jumps Y^ and Y^ with characteristics (a + 1, /3 + 1, cti, 61 — K) and (a + 1, /3, (72, b2), where cti, (T2, &i 
and &2 are defined as above. 

It is important to note that the above example have been recently used for the risk neutral stock 
price model by Eberlein and Madan |13j . 

7.3 Examples. 

Examples of Lamperti stable processes appear in the literature at least in the papers mentioned in 
the introduction ( [HI UHl [Ml HE] ) but they also appear (in a hidden way) in many other recent works. 
We will give a quick overview of some of them, not pretending to be exhaustive in this list. 
In U, we find two examples related to the factorization 

law ^fy 

e = e r„ • 

where e is an exponential variable independent of the a-stable variable r^. The first of them is related 
with the exponential functional of a killed subordinator whose Laplace exponent is given by 

^ r(aA + i) 



r(a(A-l) + l)' 

It is easy to see that it is related to the Laplace exponent <I>l of a Lamperti stable subordinator X^ 
with characteristics (a, a, a, 9), zero drift, and ct({1}) — a/r(l — a). The relationship between both 
Laplace exponents is 

0i(A)^$L(aA)+ ^^^^_^^ . 

This subordinator is also studied in Rivero |30j . where the author finds its renewal density and other 
related computations. 

The Laplace exponent of the second subordinator, here denoted by Z^, is given by 

^ r(a(A-l) + l) 
<P2(A) = A — , , — , 
1 (aA + 1) 

and can be expressed in terms of the Laplace exponent ^l,2 of a Lamperti stable subordinator AC^'^ 
with characteristics (1 — a,l,cr,9), and zero drift where (t({1}) = Q!/r(l — a). The relation between 
them is 

02(A) = Q;(f>L,2(aA). 

In both cases this allows us to compute the law of the exponential functional of aX^ and aX'^''^ 
in terms of the one oi and Z'^, respectively. 

There is another example in [l] which is related to the factorization 
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where s > a, 7s is a Gamma r.v. with parameter s and ji''^'* denotes a certain r.v which is independent 
of 7s. In this case, the killed subordinator related to the exponential functional which has the same 
moments as the 7s, can be expressed as the sum of two independent Lamperti stable procesess. In 
|30j further calculation are carried over concerning this subordinator. 

In the paper [34 in section 5.3, the authors found the Levy measure of the inverse of the local 
time at of an Ornstein Uhlenbeck process driven by a standard Brownian motion and parameter 
7 > 0. This measure is 

^3/2g7t/2 (^27)3/2e2^* 

"^^^ ~ \/2^(sinh(7t))3/2 ~ y2^(e27t - 1)3/2 

and the corresponding Laplace exponent is computed. It is related to a Lamperti stable distribution 
with characteristics (1/2, 1, \/"f/Tr)- 

This computation as well as the three former examples can be carried out by recognizing that 
behind those measures there is a related Lamperti stable distribution and applying our Theorem 1 to 
calculate the corresponding Laplace exponent. 

In the papers [10] , [54,, the main processes in study are Lamperti stable processes. All these 
papers share the property that many useful explicit calculations are be carried out. This is, we believe, 
the main advantage of this class: being at the same time a good model for many situations, allowing 
simulation of the paths as well as many explict calculation to be carried on. 
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